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ABSTRACT 

For the explosion mechanism of Type la supernovae (SNe la), different scenarios have been sug- 
gested. In these, the propagation of the burning front through the exploding white dwarf star proceeds 
in different modes, and consequently imprints of the explosion model on the nucleosynthetic yields can 
be expected. The nucleosynthetic characteristics of various explosion mechanisms is explored based 
on three two-dimensional explosion simulations representing extreme cases: a pure turbulent defla- 
gration, a delayed detonation following an approximately spherical ignition of the initial deflagration, 
and a delayed detonation arising from a highly asymmetric deflagration ignition. Apart from this 
initial condition, the deflagration stage is treated in a parameter-free approach. The detonation is 
initiated when the turbulent burning enters the distributed burning regime. This occurs at densities 
around 10 7 g cm -3 - relatively low as compared to existing nucleosynthesis studies for one-dimensional 
spherically symmetric models. The burning in these multidimensional models is different from that 
in one-dimensional simulations as the detonation wave propagates both into unburned material in 
the high density region near the center of a white dwarf and into the low density region near the 
surface. Thus, the resulting yield is a mixture of different explosive burning products, from carbon- 
burning products at low densities to complete silicon-burning products at the highest densities, as 
well as electron-capture products synthesized at the deflagration stage. Detailed calculations of the 
nucleosynthesis in all three models are presented. In contrast to the deflagration model, the delayed 
detonations produce a characteristic layered structure and the yields largely satisfy constraints from 
Galactic chemical evolution. In the asymmetric delayed detonation model, the region filled with elec- 
tron capture species (e.g., 58 Ni, 54 Fe) is within a shell, showing a large off-set, above the bulk of 56 Ni 
distribution, while species produced by the detonation are distributed more spherically. 
Subject headings: nuclear reactions, nucleosynthesis, abundances - hydrodynamics - supernovae: gen- 
eral 



1. INTRODUCTION 

There is a consensus that Type la supernovae (SNe 
la) are the outcome of a thermonuclear explosion of a 
carbon-oxygen white dwarf (WD) (e.g., Wheeler et al. 
1995; Nomoto et al. 1997; Branch 1998). For the pro- 
genitor, the Chandrasekhar-mass (Mch) WD model has 
been favored for a majority of SNe la (e.g., Hoflich & 
Khokhlov 1996; Nugent et al. 1997; Fink et al. 2007; 
Mazzali et al. 2007). 

The mass of the WD could reach Mch by several evo- 
lutionary paths, either by a mass-transfer from a binary 
giant/main-sequence companion (single degenerate sce- 
nario; e.g., Whelan & Iben 1973; Nomoto 1982) or as a 
result of merging with a binary degenerate WD compan- 
ion (double degenerate scenario; e.g., Iben & Tutukov 
1984; Webbink 1984). As the WD has accreted a suffi- 
cient amount of material, the central density of the WD 
increases and the heating rate by the carbon fusion ex- 
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ceeds the cooling rate by neutrino emission. The evolu- 
tion is then followed by a simmering phase with connec- 
tive carbon burning, lasting for about a century. At the 
end of this simmering phase, the temperature rises to the 
point where convection can no longer efficiently transport 
away the energy produced by the carbon burning. This 
is a stage where the burning becomes dynamical, initiat- 
ing a thermonuclear flame that propagates outward and 
disrupts the WD, i.e., a supernova explosion (Nomoto et 
al. 1984; Woosley & Weaver 1986). 

Once the thermonuclear flame is ignited, there are two 
possible modes of the propagation: subsonic deflagra- 
tion and supersonic detonation. A prompt ignition of 
the detonation flame is disfavored because the result- 
ing nucleosynthesis yield conflicts with Galactic chemical 
evolution (Arnett 1969) and fails to produce the strong 
intermediate-mass element features observed in SNe la. 
Thus, the explosion should start with a subsonic deflagra- 
tion flame. The deflagration stage may last until the end 
of the explosion (the deflagration model; Nomoto et al. 
1984), while it is also possible that the deflagration flame 
turns into the detonation wave [the delayed detonation 
model, or the deflagration-detonation transition (DDT) 
model; Khokhlov 1991; Yamaoka et al. 1992; Woosley & 
Weaver 1994; Iwamoto et al. 1999]. 

The supernova explosion phase has been investigated 
by "classical" one-dimensional spherically symmetric 
models; the classical deflagration W7 model of Nomoto 
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et al. (1984) has successfully explained the basic contri- 
bution of SNe la to Galactic chemical evolution, as well 
as basic features of observed spectra and light curves 
of individual SNe la of a normal class (Branch et al. 
1985). Some improvement in these observational aspects 
has been obtained by introducing a delayed detonation 
(e.g., Hoflich & Khokhlov 1996; Iwamoto'et al. 1999). 

These models, however, treat the propagation speed of 
the deflagration flame as a parameter. Moreover, the de- 
flagration flame is hydrodynamically unstable and non- 
sphericity is thus actually essential (e.g., Niemeyer et 
al. 1996). Recent investigations of the explosion models 
have intensively addressed these issues (e.g., Reinecke 
et al. 2002; Gamezo et al. 2003; Ropke & Hillebrandt 
2005; Bravo & Garci'a-Senz 2006; Ropke et al. 2006b; 
Schmidt & Niemeyer 2006). With high resolution multi- 
dimensional hydrodynamic simulations, coupled with an 
appropriate sub-grid model to capture turbulence effects 
on unresolved scales, recent studies provide essentially 
"parameter-free" simulations for the initial deflagration 
stage, where only the structure of the pre-supernova WD 
and the distribution of the initial deflagration ignition 
sparks set up the initial conditions. 

The multi-dimensional simulations have been per- 
formed with different initial conditions to sec if SNe la 
in general can be explained in a framework of a pure 
deflagration explosion, as is summarized by Ropke et al. 
(2007b). Detailed nucleosynthesis calculations have been 
performed for some deflagration models (Travaglio et al. 
2004a, 2005; Kozma et al. 2005, Ropke et al. 2006a). 
These studies indicate that the pure deflagration explo- 
sion can explain a part of SNe la, up to relatively weak 
and faint SNe la in the normal population. However, it 
has also been shown that a subsequent detonation phase 
is probably necessary to account for typical normal SNe 
la and brighter ones (Gamezo et al. 2005; Ropke et al. 
2007b). 

Delayed detonation models have also been investi- 
gated with multi-dimensional simulations (Gamezo et al. 
2005; Plewa 2007; Golombek & Niemeyer 2005; Ropke 
& Niemeyer 2007c; Bravo & Garci'a-Senz 2008). How- 
ever, compared to the deflagration models as summa- 
rized above, the investigation of multi-dimensional de- 
layed detonation models is still at the very initial stage. 
The radiation transfer based on the multi-dimensional 
models has been examined by Kasen et al. (2009). De- 
tailed nucleosynthesis studies have rarely been done (but 
see Bravo & Garcia-Scnz 2008). In this paper, we present 
results from detailed nucleosynthesis calculations, based 
on two-dimensional delayed detonation models and, for 
comparison, a pure deflagration model. 

The two delayed detonation models presented here 
can be regarded as extreme cases - not necessarily with 
respect to their 56 Ni production and brightness, but 
with respect to symmetries/asymmetries in the explosion 
phase. While in one model, the deflagration was ignited 
in an approximately spherical configuration at the cen- 
ter of the WD, the other model features an off-center 
ignition and propagation of the initial deflagration flame 
similar to the three-dimensional simulations of Ropke et 
al. (2007d). It has been suggested that the convection 
in the simmering phase may be dominated by a dipolar 
mode, and there is a good possibility that the deflagra- 
tion is initiated in an off-center way (e.g., Woosley et 



al. 2004, Kuhlen et al. 2006). In the pure-deflagration 
model, such an explosion cannot account for normal SNe 
la, because the small burning surface area should result 
in inefficient production of 56 Ni (see, e.g., Ropke et al. 
2007d). However, this can be possibly overcome in the 
delayed detonation scenario, as the detonation can po- 
tentially produce a large amount of 56 Ni. 

The paper is organized as follows. In §2, we present 
methods and models, in §3, we present our results. Dis- 
cussion of these results from a view point of the chemical 
evolution is given in §4. The paper is closed in §5 with 
concluding remarks. 

2. METHODS AND MODELS 
2.1. Explosion Models 

In this paper, we concentrate on two-dimensional mod- 
els. Although some results might be affected by the im- 
posed symmetry (e.g., Travaglio et al. 2004a; Ropke et al. 
2007b), we believe that examining 2D models is a natu- 
ral step forward. Moreover, it illustrates how nucleosyn- 
thesis in the different models proceeds, and highlights 
differences with respect to one-dimensional models. 

In particular, we focus on three models in this paper 
as follows: 

• C-DEF: A globally spherically symmetric pure- 
deflagration explosion in which the deflagration 
was ignited within the c5-shape boundary, as was 
done by Reinecke et al. (1999a). About 2 x 
10~ 2 Mq were initially incinerated to trigger the 
deflagration. This model is similar to one presented 
already in Travaglio et al. (2004a). 

• C-DDT: A delayed detonation model which fol- 
lows the 2D spherical deflagration C-DEF model. 
A prescription for the DDT is given below. 

• O-DDT: A delayed detonation model which fol- 
lows an extremely off-center deflagration. The 
deflagration is ignited by 29 bubbles distributed 
within an opening angle with respect to the z-axis 
of 45 degrees. The outermost bubble was placed 
at a distance of <~ 180 km from the center. About 
1 x 10 -5 Mq were initially incinerated to trigger 
the deflagration. The DDT is treated in the same 
way as in the C-DDT model. 

The exploding WD had a central density of 2.9 x 
10 9 gcm~ 3 . In the hydrodynamic explosion simulations, 
burning was treated in a simplified way. Only five species 
(carbon, oxygen, alpha-particles, a representative for 
iron-peak elements, and a representative for intermedi- 
ate mass elements) were followed. The deflagration and 
detonation fronts were modeled with the level-set ap- 
proach (Reinecke et al. 1999b; Golombek & Niemeyer 
2005; Roepke & Niemeyer 2007c). After passage of the 
zero-level set representing the combustion waves, the ma- 
terial was converted from the carbon/oxygen fuel mix- 
ture to an approximate nuclear ash composition. At high 
densities, NSE is reached and the ash was modeled as 
a temperature- and density-dependent mixture of iron- 
peak elements and alpha-particles. In the NSE region, 
electron captures neutronizing the ashes were followed. 
In the hydrodynamic simulations, 512 x 512 cells are used 
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TABLE 1 

Mass (M q ) and Energetic (10 51 erg) 



Form W7 C-DEF C-DDT O-DDT C-DDT/dcf. a 0-DDT/def. a 

M burn b 1.21 0.61 1.01 1.25 0.522 0.393 

E nuc c 1.78 0.91 1.46 1.80 0.767 0.522 

E K d 1.28 0.41 0.96 1.30 

E K (hyd) c 1.30 0.41 0.93 1.27 ■_ ■__ 

a At the moment when the first DDT takes place. 
^The incinerated masses. 

c The nuclear energy generation derived by the postprocess calculations. 
The asymptotic kinetic energy resulting from E nuc . The progenitor WD binding 
energy is i?bind — 0.505. 

c The asymptotic kinetic energy calculated in the hydrodynamic simulations with the 
simplified treatment of the energy generation. 



for C-DEF and C-DDT models, and 1024 x 512 cells are 
used for O-DDT model. 

In the delayed detonation models, we assumed the 
deflagration-to-detonation transition to take place once 
the flame enters the distributed burning regime (e.g., 
Ropke & Niemeyer 2007c) at a fuel density of < 1 x 
10 7 gcm -3 (hereafter /9ddt)- We emphasize that in con- 
trast to the classical ID models, the DDT is not ex- 
clusively parameterized by pddt, but the requirement 
on reaching the distributed burning implies also a tur- 
bulence criterion (see Ropke 2007a for an evaluation in 
three-dimensional models). An estimate for entering the 
distributed regime is the equality of the laminar flame 
width and the Gibson scale - the scale at which the 
turbulent velocity fluctuations equal the laminar flame 
speed. When the Gibson scale becomes smaller than the 
laminar flame width, turbulence affects (and ultimately 
destroys) the laminar flame structure. This is a prerequi- 
site for DDT, however, it is not a sufficient criterion (e.g. 
Woosley 2007; Woosley et al. 2009). The microphysics 
of DDT is not yet fully known, and therefore we apply 
the necessary condition of entering the distributed regime 
here only. We emphasize that the hydrodynamic models 
are still in a preliminary stage and used here only in or- 
der to demonstrate the nucleosynthesis associated with 
them. For more robust predictions of the 56 Ni produc- 
tion and the implied brightness of the events, more elab- 
orate hydrodynamical simulations should be used (F.K. 
Ropke et al., in prep.); the models presented here are 
understood as a case study demonstrating typical nucle- 
osynthesis for the different explosion processes. 

Results on some synthetic observables (e.g., light 
curves) derived for similar two-dimensional models can 
be found in Kasen et al. (2009), and details on hydro- 
dynamic calculations will be presented elsewhere (F.K. 
Ropke et al., in prep.). Note that the present models 
cover extreme cases in the sequence of "classical" delayed 
detonation models where the DDT takes place before the 
deflagration wave reaches the WD surface (e.g., Kasen et 
al. 2009). For the situation in which only a few initial 
bubbles are distributed within a small solid angle with 
a large off-set, the "gravitationally confined detonation" 
model (e.g. Jordan et al. 2008; Meakin et al. 2009) has 
been suggested as an alternative explosion mechanism. 
This would be a more extreme case than our O-DDT 
model. 

2.2. Nucleosynthesis 



We apply the tracer particle method to the calcula- 
tions of nucleosynthesis. The essence is to follow thermal 
histories of Lagrangian particles, which are passively ad- 
vected in hydrodynamic simulations, and then to employ 
detailed nuclear reaction network calculations to each 
particle separately. The method was first applied to core- 
collapse supernovae (Nagataki et al. 1997; Maeda et al. 
2002; Maeda & Nomoto 2003), and has become popular 
in the field thanks to its simplicity and its applicabil- 
ity to multi-dimensional problems (e.g., Travaglio et al. 
2004b). Travaglio et al. (2004a) applied the method to 
multi-dimensional, purely deflagration explosion models 
of SNe la. 

In the setup for our hydrodynamic simulations, 80 2 
tracer particles are distributed uniformly in mass co- 
ordinate, such that each particle represents the same 
mass of - 2.2 x 10" 4 M Q (= M wd /6400). The parti- 
cles are advected passively, following the velocity field 
at each time step of the Eulerian hydrodynamic simula- 
tions. The thermal history each particle experiences is 
recorded. The number of the tracer particles is sufficient 
to accurately follow the nucleosynthesis. Seitenzahl et 
al. (in preparation) find that in two-dimensional SN la 
simulations with 80 2 tracer particles all isotopes with 
abundances higher than <~ 10~ 5 arc reproduced with an 
accuracy of better than 5% (except for 20 Ne). 

The nuclear postprocessing calculations are then per- 
formed for each particle separately. To this end, we recal- 
culate the temperature from the recorded internal energy, 
rather than directly using the value obtained by the hy- 
drodynamic simulations (see Travaglio et al. 2004a). In 
deriving the temperature, the electron fraction (Y e ) is as- 
sumed to be 0.5, which introduces some errors when the 
electron captures are very active. When Tg = T/10 9 K 
> 6, we follow the abundance evolution by applying the 
Nuclear Statistical Equilibrium (NSE) abundance (i.e., 
the abundance specified by p, T, and Y c ) rather than 
fully solving the reaction network, as the NSE is reached 
in this regime of high density and temperature. In order 
to correctly follow the evolution of Y e , weak interactions 
are computed along with the NSE abundance. 

Throughout this paper, the initial C+O WD compo- 
sition is assumed as follows: X( 12 C) = 0.475, X( 16 0) 
= 0.5, and X( 22 Ne) = 0.025 in mass fractions. This is 
consistent with the compositions used in the W7 model, 
roughly corresponding to the initially solar CNO compo- 
sition: Metallicity is represented by 22 Ne, assuming that 
the CNO cycle in the H-burning has converted all heavy 
elements to 14 N, and then it is reprocessed to 22 Ne by 
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Fig. 1. — Distribution of characteristic elements/isotopes in the C-DEF model, at 10 seconds after the ignition of the deflagration wave, 
(a) Tracer particles in which the total mass fraction of stable 56 Fe, 58 Ni, and 54 Fe exceeds 0.1 (electron capture products), (b) Tracer 
particles in which the mass fraction of 56 Ni exceed 0.1, excluding particles shown in (a) (complete Si burning), (c) Tracer particles in which 
the mass fraction of 28 Si exceeds 0.1 and that of le O is below 0.1, excluding particles in (a) and (b) (O burning), (d) Tracer particles in 
which the mass fraction of le O exceeds 0.1 and that of 12 C is below 0.1, excluding particles in (a) - (c) (C burning), (e) Tracer particles 
with the mass fraction of 12 C exceeding 0.1, excluding particles in (a) - (d) (no-burning). 
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Fig. 2. — Distribution of characteristic elements/isotopes in the C-DDT model, at 10 seconds after the ignition of the deflagration wave. 
See the caption of Fig. 1. 
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Fig. 3. — Distribution of characteristic elements/isotopes in the O-DDT model, at 10 seconds after the ignition of the deflagration wave. 
See the caption of Fig. 1. 



14 N(a,7) 18 F(e+,^ ) 18 0(a, 7 ) 22 Ne in the He-burning. 

The reaction network (Thielemann et al. 1996) in- 
cludes 384 isotopes up to 98 Mo. The electron capture 
rates, which strongly affect nucleosynthesis at the begin- 
ning of the deflagration stage, are taken from Langanke 
& Martinez-Pinedo (2000) and Martinez-Pinedo et al. 
(2000). Detailed comparisons between the new rates and 
those of Fuller et al. (1982, 1985) are presented by Brach- 
witz et al. (2000) for the spherical ID deflagration model 
W7 (see also Thielemann et al. 2004). 

3. RESULTS 

3.1. Characteristic Burning Regimes 

Before presenting the result of our calculations, we 
summarize the basics of explosive burning taking place in 
the deflagration and detonation stages. Typical burning 
products can be characterized by the maximum temper- 
atures (T max ) and densities (p max ) attained by the ma- 
terial under consideration, after the passage of the ther- 
monuclear flame [see, e.g., Arnett (1996), Thielemann et 



al. (1998), and references therein, for a review of explo- 
sive nucleosynthesis] . 

The explosive burning in thermonuclear supernovae 
proceeds in different regimes, mainly characterized by the 
temperatures reached in the nuclear ashes (Thielemann 
et al. 1986). In deflagrations, this directly translates into 
characteristic fuel densities ahead of the flame. At the 
first stage of the deflagration, the density is higher than 



Temperatures rise to T maXi g = T max /10 9 K 



10 8 g cm 3 . 

> 6. At this temperature, NSE applies, and thus the final 
composition is determined by freezeout processes (repre- 
sented by T max and /9 max ) and the efficiency of electron 
capture reactions. Because of the high density, electron 
capture reactions are important. Dominant species in 
this "complete silicon burning with electron capture re- 
gion" are stable 56 Fe, 54 Fe, 58 Ni, and radioactive 56 Ni 
(which decays into 56 Fe). In the case of even stronger 
electron capture reactions, the main products are 50 Ti, 
54 Cr, and 58 Fe (e.g., Thielemann et al. 2004). 
Following the expansion of the WD, the deflagration 
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Fig. 4. — Temporal evolution of temperature distribution (at 0.7, 0.9, 1.1, 1.3, 1.5, 1.7 seconds after the ignition of the deflagration flame, 
from left to right), in the O-DDT model. The color coordinate corresponds to Tg = 8 (red) and (blue). The complete silicon burning 
takes place at the temperature above Tg = 5 (yellow). 
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Fig. 5. — Temporal evolution of distribution of NSE elements (at 0.7, 0.9, 1.1, 1.3, 1.5, 1.7 seconds after the ignition of the deflagration 
flame, from left to right), in the O-DDT model. Shown here is the mass fraction, from (blue) to 1 (red). 



proceeds progressively at the lower density. At T ma x,9 ~ 
5 - 6 and /w ~ 5 x 10 7 - 10 8 g cm- 3 , NSE still' ap- 
plies. The electron captures are no longer important, and 
thus the initial Y c is virtually preserved. The dominant 
species are 56 Ni and 58 Ni. 

In contrast to spherically-symmetric delayed detona- 
tion models, the detonation in our two-dimensional se- 
tups is triggered at certain spots at the flame front, but 
at other locations, deflagration burning can still proceed 
for some while. It is characterized by successively lower 
temperature and density. At T max 9 ^3 — 5, oxygen 
burning or incomplete silicon burning is a result, leaving 
mainly intermediate-mass-elements (IME) such as 28 Si 
and 32 S. At T max> g ^2 — 3, carbon burning or neon 
burning is a result, characterized by abundant 16 and 
24 Mg in consumption of 12 C. At even lower temperatures 
(T maXi g < 2), virtually no major thermonuclear reactions 
take place. 

The detonation wave can also be characterized by the 
same tendency, although the temperature reached in the 
ashes is not a unique function of the fuel density any- 
more but also depends on the shock strength, the effect of 
which needs further investigation in hydrodynamic sim- 
ulations. In the delayed detonation framework, the det- 
onation is ignited after significant expansion of the WD, 
and thus the electron captures are almost always unim- 
portant (§3.4). If there is still a high density region left 
after the deflagration, the detonation wave can convert 
the material chiefly into 56 Ni by complete silicon burning. 
As the density decreases, the detonation nucleosynthesis 



is characterized by oxygen burning, carbon burning, and 
eventually burning ceases. Compared to the correspond- 
ing deflagration case at a similar burning stage, these 
stages are usually encountered at lower fuel densities, as 
the detonation compresses the burning material. 

3.2. Structure of the ejecta 

Table 1 shows the global features of our calculations. 
In the C-DEF model, more than half of the WD mate- 
rial is left unburned. The final kinetic energy, ~4x 10 50 
erg, is significantly lower than that inferred for normal 
SNe la. In the C-DDT model, ~ 0.5M Q of the C/O 
material is incinerated by the detonation wave, after the 
deflagration already has burned ~ O.5M . The inciner- 
ated mass and nuclear energy release in the deflagration 
stage of the O-DDT model are smaller than the C-DDT 
model. The detonation wave, however, burns a larger 
amount of material and produces a larger amount of nu- 
clear energy in the O-DDT model than in the C-DDT 
model. The final incinerated mass and the kinetic en- 
ergy in the O-DDT model are comparable to those in 
the W7 model. Hereafter, we discuss how these results 
can be understood in terms of the flame propagation and 
the DDT. Note that the nuclear energy release in the hy- 
drodynimic simulations (with a simplified treatment for 
the nuclear reactions) and that in the detailed reaction 
network calculations are consistent within 5%. 

Figures 1-3 show the distribution of selected species 
at 10 seconds after the ignition of the deflagration flame. 
At this time, burning has ceased and the SN material is 



6 



E 1° 



- «— 1 

{a) C-DEF 




_>_i 


■ 

i . . 



E 10 




E 10 




T .„ IK] 



TJ(10 K} 



Fig. 6. — Maximum density p max and temperature T u 
flame. The solid line is for the W7 model. 



for individual tracer particles (dots), after the passage of the thermonuclear 



already almost in a homologous expansion. The C-DEF 
model (Fig. 1) shows large-scale mixing, and thus it does 
not possess a clearly layered structure as predicted in 
exactly spherical ID models (e.g., W7). A large amount 
of unburned carbon and oxygen are left, being mixed 
down to the central region. 

In contrast, the detonation wave produces a more 
layered structure, since this supersonically propagating 
wave is unaffected by hydrodynamic instabilities (at least 
on the large scales resolved in the models). In the C-DDT 
model, the deflagration burns out much of the center of 
the WD before the detonation is triggered. Therefore, 
the detonation mainly burns material towards the stel- 
lar surface, but it also propagates down the fuel funnels 
between the fingers of the deflagration ash (Fig. 1). In 
the C-DDT model, the whole WD, including the central 
region, experiences the strong expansion in the vigorous 
deflagration stage. Therefore, the unburned material left 
after the deflagration is processed in the detonation stage 
mainly by carbon and oxygen burning, producing O and 
IME, but virtually no Fe-peak elements (Fig. 2). In this 
particular model, Fe-peak elements are thus produced 
mostly in the initial deflagration stage. 

In the O-DDT model the burning proceeds in an as- 
pherical manner. Figures 4 and 5 show the temporal 
evolution in the O-DDT model. The deflagration frame, 
ignited off-center, floats outward, and spreads laterally 
(Ropke et al. 2007d). It creates a large blob of ashes in 
the upper hemisphere of the WD. Due to the lateral ex- 
pansion in the outer layers of the star, the neutron-rich 
Fe-peak elements produced by electron captures end up 
in a characteristic off-center shell- like layer (Fig. 3). The 
detonation triggers on top of this blob at ~ 1 second. It 
cannot cross ash regions (Maier & Niemeyer 2006) and 
has to burn around the ash blob in order to reach the 
central parts of the WD. Consequently, it initially prop- 
agates only outward (at ~ 1.1 second). 

At about 1.5 sec after the deflagration ignition, the 
detonation wave has burned around the ash blob and 
propagates inward. Since the off-center deflagration did 
not release much energy, the central parts of the WD are 
still dense and contain mostly unburned material. This 
material is converted to Fe-peak (predominantly 56 Ni) 
by the detonation wave. 

3.3. Characteristic thermal properties 

To understand the characteristic structure of the ejecta 
described in §3.2, we examine details of the thermal his- 
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Fig. 7. — Average peak density of tracer particles as a function of 
peak temperature (with the temperature bin taken to be AT max = 
2 X 10 s K). Models are W7 (blue thin-solid), C-DEF (green thin- 
dashed), C-DDT (magenta thick-dashed), and O-DDT (red thick- 
solid). 

tory of the tracer particles. Figure 6 shows the max- 
imum temperature (T max ) and density (/) max ) for each 
tracer particle (§3.1) after flame passage. For compari- 
son, the same values for the classical W7 model are also 
shown. In Figure 7, average behaviors are plotted, which 
were obtained by taking the mass-average as a function 

of ^max- 

Figures 6 and 7 show that the distribution of (T max , 
Pmax) is similar for all the models at T maXj 9 > 3 (corre- 
sponding to pmax ~ Pddt). This reflects the fact that 
the local properties of the deflagration flame are basi- 
cally independent from the geometry. Interestingly, de- 
spite the different prescription of the deflagration flame 
propagation and different expansion time scale (see e.g., 
Travaglio et al. 2004a) , we note a similarity with the out- 
come of the W7 in this respect. 

At T maXj g < 3, the models differ. Note that the be- 
havior of W7 in this temperature regime shown in Fig- 
ures 6 and 7 is not real; the deflagration is turned off 
at Tmax, 9 ~ 2, where the flame is already near the sur- 
face in the W7 model. Our deflagration model extends 
to lower T max and /o max . In the two delayed detonation 
models, this temperature range corresponds to the deto- 
nation phase, and thus the behavior is expected to devi- 
ate from the deflagration models, and indeed it does. The 
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Fig. 8. — The masses (Mg) within each peak temperature range 
(AT = 2 X 10 8 K). Models are W7 (blue thin-solid), C-DEF (green 
thin-dashed), C-DDT (magenta thick-dashed), and O-DDT (red 
thick-solid). 

bution; the compression and the maximum temperature 
are roughly linear functions of the unburned density. 

Figure 6 shows a significant dispersion in the T max — 
Pmax distribution around the mean value, reaching nearly 
one order of magnitude in p max for given T max , at 
T maXj g < 5. This is a unique behavior in the multi-D 
models, potentially leading to a variety of burning prod- 
ucts otherwise not expected in strictly spherically sym- 
metric models like W7. 

The different hydrodynamic properties lead to differ- 
ent amounts of material as a function of T max (or p max ), 
and this is the main reason why different nucleosynthesis 
features appear. Figure 8 shows the mass of the material 
as a function of T max . It is seen that the deflagration 
is efficient in W7, and thus the W7 model has a much 
larger amount of the mass with T max g > 6 than the 2D 
models. At T maXi g < 3 — 4, there appears a peak in the 
C-DDT model; this is due to the detonation wave burn- 
ing the fuel material left behind the deflagration stage. 
The O-DDT model shows a characteristic distribution; 
two peaks at T maXi g ~ 3 — 4 and at T max ,g ^5 — 6. The 
former is the detonation wave propagating outward, as 
also seen in the C-DDT model. The latter one at higher 
Tmax is due to the detonation wave propagating inward, 
near the central region where the fuel still has high densi- 
ties. The gap between the two peaks corresponds to the 
(off-center) shell structure of the deflagration products 
(i.e., neutron-rich Fe-peak elements; Fig. 3). 

3.4. Deflagration and Electron Captures 

In the densest region in the deflagration stage, electron 
capture reactions proceed to synthesize neutron-rich Fe- 
peak elements under NSE conditions. Figure 9 shows 
the electron fraction (Y e ) as a function of /0 max , for the 
O-DDT model. The same figure for the other 2D models 
(not shown) is similar to the O-DDT model. Figure 9 
shows that the electron capture becomes important at 
pmax >5x 10 8 g cm -3 , and Y e as low as 0.463 is realized 
in the highest density region. Thus, the main product of 
the highest density region is stable 56 Fe (Y e ~ 0.464), fol- 
lowed by 54 Fe (Y e ~ 0.481) and 58 Ni (Y e ~ 0.483) in the 
lower density region. Since Y c > 0.46 , 50 Ti (Y c ~ 0.440), 
54 Cr (Y c ~ 0.444) and 58 Fe (Y c ~ 0.448) are not pro- 



Fig. 9. — The final values of Y e as a function of the peak density, 
for the O-DDT model. 

duced abundantly in the present models. The behavior 
in Y c and the resulting electron capture reactions are 
largely consistent with the 2D deflagration model pre- 
sented in Travaglio et al. (2004a). The electron cap- 
tures, however, are less efficient in the 2D models than 
in the 3D deflagration model (Travaglio et al. 2004a) and 
in the classical ID models (Thielemann et al. 2004; see 
also §3.6). This is likely due to the weaker deflagration 
in the present models than in the others (ID and 3D), 
although an uncertainty is involved in the treatment of 
the temperature evolution in the NSE regime (§2.2). 

Figure 10 shows the comparison between Y e before and 
after the DDT. In the C-DDT model, Y c is not affected 
by the detonation wave as is consistent with previous ID 
models. On the other hand, in the O-DDT model the 
detonation slightly affects Y c especially of the material 
which was hardly processed by the deflagration (~ Y c > 
0.496). The change in Y c introduced by the detonation 
is at most AY C ~ 0.02, much smaller than the change 
in the deflagration stage at > 10 8 — 10 9 g cm -3 . This is 
a result of electron capture reactions at the detonation 
wave propagating near the WD center, since the central 
region is still at high density (e.g., Meakin et al. 2009). 
The change is not as large as that in the GCD models 
of Jordan et al. (2008) and Meakin et al. (2009), which 
resulted in AY e ~ 0.05. This is because the distribution 
of the initial bubbles in our simulations is not as extreme 
as theirs (see §2.1) and the WD pre-expansion before the 
onset of the detonation was stronger in our model. 

3.5. Distribution of Nucleosynthesis Products 

Figure 11 shows the radial abundance distribution, for 
which the material is angle-averaged within each veloc- 
ity bin. Note that our models are followed up to ~ 10 
seconds, thus the velocity is proportional to the radial 
distance, and the distribution in velocity space corre- 
sponds in good approximation to that in spatial space. 
The O-DDT model has a very aspherical ejecta struc- 
ture. Therefore, the angle-dependent radial abundance 
distribution in this model is given in Figure 12. 

In the C-DEF model, different burning regions are 
macroscopically well-mixed. As a consequence, it forms 
a broad region in which the electron capture products 
( 56 Fe, 54 Fe, and 58 Ni), complete-silicon burning products 




Fig. 10. — The final values of Y e (i.e., after the detonation) as compared to those just before the first DDT transition (i.e., before the 
detonation), for (a) C-DDT and (b) O-DDT models (shown by dots). 



( 56 Ni), intermediate mass elements ( 32 S, 28 Si, 24 Mg), and 
unburned elements ( 16 and 12 C) coexist (when averag- 
ing in the zenith angle). This region is surrounded by 
unburned material at radii not reached by the flame. An 
enhancement of C and O near the center due to the large- 
scale deflagration mixing is visible. These are typical of 
2D deflagration models; this smoothing effect tends to 
be suppressed in 3D deflagration models (Ropke et al. 
2007b). 

In the C-DDT model, the inner region below ~ 8, 000 
km s~ x has a similar structure to that in the deflagration 
model. The difference is seen in the innermost region, 
where unburned C and O are burned into intermediate 
mass elements by the detonation wave. In this model, 
C is almost completely burned, but there still remains 
O with a mass fraction of ~ 0.1, since the detonation in 
this model leads to carbon or oxygen burning. 

Surrounding the deflagration region is the oxygen 
burning layer, where 28 Si, 32 S, and 24 Mg are the main 
products. This is surrounded by a mixture of carbon 
burning products and unburned material. Note that the 
carbon burning region and the unburned region are mi- 
croscopically separated. As shown in Figure 2, the un- 
burned region (with almost initial C+O composition) is 
not fully homogeneously distributed, and there are "fin- 
gers" of carbon burning regions (where oxygen and IMEs 
are produced in consumption of carbon). 

The radial composition distribution in the O-DDT 
model follows that of the C-DDT model. A striking 
difference, however, is that the region of deflagration 
ashes is not in the center but confined within V ~ 
5,000 — 10,000 km s _1 . The innermost region is con- 
verted to 56 Ni, by the complete silicon burning taking 
place at the inward detonation wave. 

3.6. Total Yields 

Figure 13 shows the nucleosynthesis yield for each 
model, integrated over the whole ejecta. The ratios of 
the masses of elements to those of the W7 model are 
plotted in Figure 14. Tables 2 and 3 list the masses of 
the elements and isotopes (after radioactive decays) while 
major long-lived radioactive isotopes are given in Table 
4. For the C-DDT and O-DDT models, the values at the 
time of the first DDT (~ 1.15 sec for C-DDT and - 1.05 
sec for O-DDT) are also shown in Tables 2 and 3. 



For comparison, we present the yields of the W7 model, 
calculated by ourselves using the thermal history of the 
original model. Details are different from Iwamoto et 
al. (1999), because of the updated electron capture rates 
(see Brachwitz et al. 2000). 

The result for the C-DEF model is consistent with the 
similar 2D model in the previous study (Travaglio et al. 
2004a). First of all, the mass of 56 Ni is only - 0.25M Q , 
smaller than in the W7 model (~ 0.64A/©). A large 
amount of C and O are left unburned. The final C/O ra- 
tio is larger than in the W7 model, because of the weak 
C-burning in our model resulting in almost the original, 
unburned C/O ratio near the surface (Fig. 11). IMEs 
are underproduced as compared to the W7 model. Note, 
however, that generally full three-dimensional models al- 
leviate this problem (e.g., Travaglio et al. 2004a; Ropke 
et al. 2007b). 

In the C-DDT model, the density of the WD is already 
low when the DDT takes place, and the temperature of 
material after the passage of the detonation wave does 
not reach ~ 5 x 10 9 K (Fig. 8). The detonation therefore 
produces at most IMEs; Fe-peak elements are produced 
in the initial deflagration stage, and the mass of 56 Ni is 
almost the same as in the deflagration model. Conse- 
quently, the abundance pattern of the Fe-peak elements 
(including 56 Ni) is very similar to that of the C-DEF 
model (Figs. 13 and 14). IMEs are produced much more 
abundantly than in the W7 model. Because of the carbon 
burning, the C/O ratio is much smaller than for W7. 

Somewhat surprisingly, the O-DDT model predicts the 
abundance pattern similar to W7 (Fig. 13), with most 
of elements consistent with W7 within a factor of two 
(Fig. 14). The mass of 56 Ni (~ 0.54Af Q ) is compara- 
ble with the W7 model (~ O.64M ) (but see §4). IMEs 
tend to be overproduced by a factor of ~ 2, because of 
the carbon and oxygen burning in the detonation phase. 
Carbon is almost entirely consumed by the carbon burn- 
ing- 

A comparison between the W7 model and our 2D mod- 
els shows that very strong electron capture reactions 
are not efficient in the 2D models (§3.4). 50 Ti, 54 Cr, 
and 58 Fe, produced at the highest density resulting in 
Y c < 0.46, are not abundantly produced in the 2D mod- 
els. 



9 





Fig. 11. — Mass fractions of selected isotopes as a function of the expansion velocity, averaged in the velocity bin of AV. (a, d) C-DEF 
model, AV = 250 km s" 1 , (b, e) C-DDT, AV = 500 km s~\ and (c, f) O-DDT, AV = 500 km s" 1 . 
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Fig. 12. — Angle-dependent mass fractions of selected isotopes in O-DDT, as a function of the expansion velocity, averaged in the velocity 
bin of AV = 500 km s _1 . (a, d) averaged within 0-37 deg, (b, e) 78 - 90 deg, and (c, f) 143 - 180 deg. 
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Fig. 13. — Isotopic yields integrated within the whole ejecta (after radioactive decays). The mass fraction is shown relative to the solar 
value, and normalized by 56 Fe. 
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Fig. 14. — Ratios of integrated element yields to the W7 model yields. 
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TABLE 2 

Masses of elements (Mq) 
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a At the moment when the first DDT takes place. 



3.7. Radioactive Isotopes 

The amount of some radioactive isotopes in our deto- 
nation models is quite different from W7 and even with 
canonical, spherical delayed detonation models. Com- 
parison between our results (Tab. 4) and Table 4 of 
Iwamoto et al. (1999) shows that radioactive isotopes 
lighter than Fe-peak tend to be more abundant, with 
the tendency of larger amounts of these isotopes for 
smaller density at the detonation. On the other hand, the 
amount of some Fe-peak radioactive isotopes can change 
dramatically from the classical ID delayed detonation 
model, due to updated electron capture rates in our cal- 
culations (§2.2). 

4. DISCUSSION: IMPLICATIONS FOR CHEMICAL 
EVOLUTION 

In this paper, we present our first results for nucle- 
osynthesis in 2D delayed detonation models. The models 
investigated in this paper represent extreme cases 6 ; the 
deflagration is initiated either at the center (C-DDT) or 
by off-center bubbles confined in the narrow angle with 
respect to the z-axis (O-DDT). The important variation, 
we have not examined in this paper, is the deflagration 
bubbles distributed more or less spherically, but at a dis- 
tance from the center (e.g., Kasen et al. 2009). In this 
case, we expect that the propagation of the burning will 
produce an ejecta structure that falls in between the O- 
DDT and the C-DDT models. Some material in the cen- 
tral region is left unexpanded in the deflagration stage, 
and the detonation can burn the material there to the 
Fe-peak elements. Thus, we expect that nucleosynthe- 
sis features are similar to the O-DDT model, except for 
the strong angle-dependence seen in the O-DDT model. 

6 Note that the O-DDT model is still less extreme than the CCD 
model (Jordan et al. 2008; Mcakin et al. 2009) in the distribution 
of the initial bubbles. 



Depending on the initial placement of the deflagration 
bubbles, the deflagration ashes may stay more confined 
to the center than in the O-DDT model. 

Also, we should note that the treatment of the DDT in 
the hydrodynamic calculation is still preliminary (§2.1). 
According to simulations with several different prescrip- 
tions, especially with different values for pddt, we be- 
lieve that the prescription used in this paper results 
in a relatively weak detonation (see also Gamezo et al. 
2005; Bravo & Garcfa-Scnz 2008). Therefore, although 
the models presented here are extreme with respect to 
the ejecta structure, they do not necessarily reflect the 
range of 56 Ni production possible within the framework 
of the delayed detonation model. The results of an ex- 
tended survey of delayed detonation models will be pre- 
sented elsewhere (F. Ropke et al., in prep). Despite these 
caveats, the general features found in this paper are not 
expected to be sensitive to such details. 

The nucleosynthesis features in the W7 model have 
been shown to be roughly consistent with the Galactic 
chemical evolution (e.g., Iwamoto et al. 1999; Goswami & 
Prantzos 2000). One problem is that the ratio 58 Ni/ 56 Fe 
is too large in the W7 to compensate the over-solar pro- 
duction of 58 Ni in core-collapse SNe. The problem is ba- 
sically solved by introducing a delayed detonation. Af- 
fecting primarily lower density material, it produces a 
large amount of 56 Ni but virtually no 58 Ni due to the in- 
efficiency of electron captures here. We see that the ratio 
58 Ni/ 56 Fe is slightly decreased in the O-DDT model com- 
pared to the W7, marginally satisfying the constraint 
from the chemical evolution [( 58 Ni/ 56 Fe)/( 58 Ni/ 56 Fe) 
~ 3]. The situation should become better if the deto- 
nation burns more material to NSE than in the present 
prescription. We should also note one limitation in the 
present analysis, i.e., the treatment of the temperature in 
the NSE regime. The temperature is extracted assuming 
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TABLE 3 
Masses of isotopes (M©) 
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TABLE 4 

Masses of radioactive species (M©) before 

RADIOACTIVE DECAYS. 



Isotopes 
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Y e = 0.5, which should introduce some errors when the 
electron capture reactions are active. This could to some 
extent affect the nucleosynthesis in the initial deflagra- 
tion phase (§3.4), and therefore the ratio 58 Ni/ 56 Fe. 

The intermediate mass elements can also be used 
to constrain SN la models. Assuming a typical frac- 
tion of ~ 20% for the frequency of SNe la as com- 
pared to core-collapse SNe, then (Si/Fe)/(Si/Fe) < 0.5 
is required to compensate the over-solar production in 
core-collapse SNe. The O-DDT model has the ratio 
(Si/Fe)/(Si/Fe)© ~ 1; thus, the model can explain at 
most a half of SNe la. The C-DDT is not at all fa- 
vored, since (Si/Fe)/(Si/Fe) Q ~ 2, i.e., it produces too 
much IMEs relative to Fe. The deflagration models (W7 
and C-DEF) do not have the problem in the Si/Fe ratio 
[(Si/Fe)/(Si/Fe) ]Q < 0.5]. In the O-DDT model, all the 
elements heavier than N are consistent with W7 within 
a factor of two. 

In reality, the solar abundances are the superposition 
of contributions from SNe at various metallicities, and we 
expect that the " average" SNe la to occur at a metallicity 
less than solar. We indeed expect some improvement if 
we consider a lower metallicity (see e.g., Travaglio et al. 
2005): a decreased metallicity should result in a smaller 
58 Ni/ 56 Fe ratio in the delayed detonation models. Also, 
the smaller metallicity is expected to lead to the smaller 
amount of 54 Fe. Mg and Al are also affected: If the 
metallicity is 10% of the solar value, then the amount of 
24 Mg could increase by ~ 50%, and that of 27 A1 could 
decrease by a factor of a few. These changes do not 
conflict with the Galactic chemical evolution, since the 
ratios (Mg, Al) /Fe are much smaller than the solar values 
in the O-DDT models. 

Summarizing, the hypothesis that about half of SNe 
la are represented by the extreme O-DDT model, is not 
rejected by the chemical evolution arguments. Note that 
this should also apply to globally symmetric, but off- 
center, delayed detonation models to some extent, as 
these models should share the basic feature that the cen- 



tral region is burned to Fe-peak and 56 Ni by the detona- 
tion (§4). Therefore, multi-D delayed detonation models 
can potentially account for a majority of SNe la. 

5. CONCLUDING REMARKS 

We have presented results of the detailed nucleosynthe- 
sis calculations for 2D delayed detonation models, focus- 
ing on an extremely off-center model (O-DDT model). 
Features are different from classical ID models. Unlike a 
globally symmetric delayed detonation model following 
the central ignition of the deflagration flame (C-DDT 
model), the detonation wave proceeds both in the high 
density region near the center of a white dwarf and in the 
low density region near the surface. Thus, the resulting 
yield is a mixture of different explosive burning prod- 
ucts, from carbon-burning at low densities to complete 
silicon-burning at the highest densities, as well as the 
electron-capture products from the deflagration stage. 

The evolution of the deflagration flame can be different 
in 2D and in 3D simulations (e.g., Ropke et al. 2007b). 
We believe that the global feature found in this study, 
i.e., the detonation wave propagating into the innermost 
region, would not be changed substantially in 3D sim- 
ulations. However, some details would be affected: For 
example, 3D simulations tend to result in stronger de- 
flagration than in 2D simulations for similar initial con- 
ditions (e.g., Travaglio et al. 2004a). As the stronger 
deflagration is followed by the weaker detonation in our 
DDT models, this would result in a shift of the overall en- 
ergetic and the nucleosynthesis production as compared 
to the 2D models. On the other hand, although in our 
2D simulations the detonation wave has to burn around 
the deflagration ash blob before reaching to the center, 
the deflagration ash blob may have holes, depending on 
the ignition geometry through which the detonation wave 
can directly penetrate into the center in 3D simulations. 
Thus, future 3D simulations are important to obtain the 
robust model predictions as a function of the model input 
(e.g., distribution of the deflagration ignition sparks). 

The yield of the O-DDT model largely satisfies con- 
straints from the Galactic chemical evolution despite the 
low DDT density as compared to ID delayed detonation 
models. This is a result of qualitatively different behavior 
of the detonation propagation in ID and multi-D mod- 
els; the outward propagation in former (which is also the 
case in the C-DDT model), while the inward propagation 
in the latter. 

The O-DDT model could thus account for a fraction of 
SNe la (especially bright ones). As less-extreme (more 
spherical) models are expected to share the common 
properties in the integrated yield, the multi-D delayed 
detonation model could potentially account for a main 
population of SNe la. 

The delayed detonation models also provide a charac- 
teristic layered structure, unlike multi-dimensional, espe- 
cially 2D, deflagration models 7 . In the O-DDT model, 
the region filled by electron capture species (e.g., 58 Ni, 
54 Fe) shows a large off-set, and the region is within a shell 
above the bulk of the 56 Ni distribution near the center. 
These can be directly tested by observations. 

We note that the distribution of the nucleosynthsis 

7 Note that the 3D deflagration models can also potentially pro- 
duce the layered structure (Ropke et al. 2007b). 
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products in our 2D DDT models is somewhat different 
from the result of Bravo & Garcia-Senz (2008). Their 
3D DDT models lack a clear abundance stratification, 
and they are characterized by large mass fractions of Fe- 
peak elements near the surface regions. In contrast, the 
3D DDT models of Ropke & Niemeyer (2007c), do pro- 
duce a layered structure with the surface dominated by 
IMEs - similar to our present 2D models (§3.5). It seems 
that (1) the deflagration is stronger in Bravo & Garcia- 
Sent (2008) than in our 2D models and in 3D models of 
Ropke & Niemeyer (2007c), and (2) the detonation wave 
is mainly propagating inward in Bravo & Garcia-Senz 
(2008) although it is propagating both inward (producing 
Fe-peak elements) and outward (producing IME) in our 
simulations. As a result, the amount of unburned mate- 
rial is smaller in our models. The cause of the different 
flame propagation is not clear, but likely due to different 
treatment of thermonuclear flames. The overall abun- 
dance distribution of the 3D DDT models of Gamezo 
et al. (2005), on the other hand, is similar to our 2D 
models, producing the stratified configuration. In their 
simulations, they initiated the detonation from the cen- 
ter at relatively high DDT density, and the detonation 
wave propagates outward, producing the layered compo- 
sition structure as the temperature drops following the 
detonation propagation. This is, indeed, quite different 
from our models, in which the detonation is initiated at 
relatively low DDT density but propagates inward to the 
high density central region. 

The observational consequences from similar models 
have been discussed by Kasen ct al. (2009) for the 
early photospheric phase (see also Hillebrandt et al. 2007; 
Sim et al. 2007). Here, we summarize some additional, 
expected observational characteristics, especially in the 
late-time nebular phase (taken after a few hundred days) . 
We emphasize that the late-time spectroscopy is cur- 
rently the most effective way to hunt for the signature 
of the DDT model in the innermost region. See also 
Maeda et al. (2010) who discussed the following points 
in details. 

• Abundance Stratification and carbon near 
the surface: Overall, the stratified composi- 
tion structure obtained in our 2D models is con- 
sistent with the result of the "abundance tomog- 
raphy" (Stehle et al. 2005; Mazzali et al. 2008). 
Spectrum synthesis models as compared to the ob- 
served early-phase photospheric spectra show that 
the mass fraction of carbon at ~ 10, 000 — 14, 000 
km s _1 should be smaller than 0.01 (e.g., Branch et 
al. 2003; Thomas et al. 2007; Tanaka et al. 2008). 
This constraint is satisfied by our DDT models at 
10,000 km s" 1 , although it is marginal at 14, 000 
km s _1 . The latter could, however, be improved 
by changing the DDT criterion such that the tran- 
sition on average proceeds at higher densities (e.g., 
Iwamoto et al. 1999). Another interesting observa- 
tional target is the inhomogeneous distribution of 
the unburned C+O pockets near the surface both 
in C-DDT and O-DDT models (§3.5), although the 
total amount of such unburned material is small in 
the DDT models. Detectable carbon absorption 
lines may appear only when the observed linc-of- 
sight intersects a large number of such unburned 



pockets, which may explain the low frequency of 
SNe la showing the C absorption lines detected 
(e.g., Tanaka et al. 2006). 

• [O I] AA6300,6363: The main problem in the 
pure deflagration model is existence of unburned 
carbon and oxygen mixed down to the central re- 
gion. This should produce a strong [O I] AA6300, 
6363 doublet in late-time spectra, although no such 
signature has been detected in the observations 
(Kozma et al. 2005). This tension is alleviated in 
3D deflagration models, with the mass fraction of 
unburned elements going down to 10 % (Ropke et 
al. 2007b) or even to smaller values (e.g., Travaglio 
et al. 2004a). The DDT model does not have this 
problem, as unburned material near the center is 
burned in the subsequent detonation. A detailed 
study of late-time nebular spectra should provide 
us with information on this issue. 

• Line shifts: Profiles of nebular emission lines 
can be used to effectively trace the distribution of 
the burning products, as is proven to be efficient 
for core-collapse SNe (Maeda et al. 2002; Mazzali 
et al. 2005; Maeda et al. 2008; Modjaz et al. 2008; 
Taubenberger et al. 2009). The O-DDT model pre- 
dicts that (1) the distribution of stable Fe-peak el- 
ements is off-set following the asymmetric deflagra- 
tion flame propagation, while the detonation prod- 
ucts (e.g., a large fraction of 56 Ni) are distributed 
more or less in a spherical manner (Figs. 5 and 
12). Recently, Maeda et al. (2010) found that the 
expected variation of the line wavelength is seen 
in nebular spectra of SNe la, indicating that the 
above configuration can be relatively common in 
SNe la. 

• Line profiles: Detailed profiles of the nebular 
emission lines can be used to infer the distribution 
of emitting ions. It has been suggested to use Near- 
Infrared (NIR) spectroscopy to hunt for the geom- 
etry, as NIR [Fe II] lines are not severely blended 
(e.g., Hoflich et al. 2004; Motohara et al. 2006). 
This can also be done using several IR lines like 
[Co III] 11.88/im (Gerardy et al. 2007). The best 
data to date have been obtained for SN 2003hv, 
showing a hole in the distribution of 56 Ni 8 . Be- 
cause the mixing is introduced by the deflagration, 
such a hole is difficult to understand in the present 
models (and most of SN la explosion models; but 
see Meakin et al. 2009). This issue remains unre- 
solved, and need further study preferentially in 3D 
simulation. 
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